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It is currently believed that the decay of the thermoremanent magnetization in spin glasses is 
composed of two terms; The "stationary" term which does not depend on the sample history and 
dominates the short time decay (< Is) and a long time aging term which depends on the samples 
history. The sample history includes both the thermal history and the time spent in a magnetic 
field (waiting time) before that field is removed. We report finding a third component of the 
decay at times much larger than the waiting time. This decay is independent of the waiting time, 
but part of the same mechanism that produces aging. We explain this decay in terms of the 
initial state distribution that is developed during the cooling process. This explanation is tested by 
performing Iso-thermoremanent magnetization decay measurements. Follow up calculations using 
the Spin Glass phase space barrier model suggest that the long time decay is logarithmic and that 
the maximum aging time is approximately three thousand years, in this sample, at a measuring 
temperature of .83 Tg. 

PACS numbers: 75.50 Lk 



Understanding the effects of disorder on states of mat- 
ter has been the subject of intense study for almost 100 
years. In the early 1970s the realization that highly dis- 
ordered magnetic materials with random interactions un- 
dergo an apparent phase transition0,|3|, sparked consid- 
erable interest and excitement in the condensed matter 
community. This discovery of the so-called spin glass 
phase was marked by a large experimental effort by many 
researchers in an attempt to try to understand the phys- 
ical parameters of the underlying phase spaceQ- One 
very interesting property of the spin glass phase was the 
observation that there exists a large distribution of relax- 
ation times extending from atomic fluctuation time scales 
through what are believed to be geological time scales. 

In 1983 it was shown that materials in the spin glass 
phase have memory effects on time scales ranging from 
milhseconds, to seconds, to days and beyond0|. To date 
much information has been obtained about the structure 
of the spin glass phase space from the study of aging ef- 
fects. The classic measurement of these memory effects 
is the Thermoremnant (TRM) or complementary Zero 
Field Cooled (ZFC) magnetization decay measurements. 
In the TRM experiments the sample is cooled, in a small 
constant magnetic field, through its transition temper- 
ature, to a measuring temperature T^- After waiting a 
time im, in the magnetic field, the field is rapidly removed 
and the consequent magnetization decay of the sample 
measured. The decay is strongly dependent on the wait- 
ing time itw)- This is called aging. It has recently been 
observedpj that the magnetization decays scale with tw 
for a spin glass cooled, sufficiently rapidly, to a particular 



measuring temperature. 

While aging has been the primary focus of TRM mea- 
surements, it is not the only contribution to the decay. 
Upon cooling the spin glass sample from a high temper- 
ature, in a magnetic field, the magnetization generally 
follows a Curie or Curie Weiss like behavior. Cooling 
the spin glass further, through its transition temperature, 
the magnetization appears to remain approximately con- 
stant, at a value Mjc- This is an indication that the spins 
freeze in approximately random directions with a net bias 
due to the field. Upon removal of the field there is a rapid 
decrease in the sample magnetization, from a value of 
Mfc^ before entering the aging regime. Various measure- 
ments including AC susceptibility [1, 0, muon decay 
neutron spin echoQ and NMRjl^l each provide some ev- 
idence that this rapid decay is power law in nature. The 
power law term is independent of the waiting time and 
has been called the stationary term (i.e. stationary with 
respect to the waiting time). 

Current belief 0,0,0] is that the full magnetization 
decay can be described by a power law plus an aging term 
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In this study we show that there exists another part of 
the decay which is independent of the waiting time and 
occurs in the time regime t » tw This is a new region 
of the decay and is not related to previous discussions 
of t >> ito0, 0|. Analysis of this long time regime 
suggests that it is independent of the waiting time and is 
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FIG. 1: A(a) Decay curves for effective waiting times 
of: 7,17.5, 25, 27, 35, 44 seconds. In the inset, the 
S(t)'s for each curve used to determine the effective wait- 
ing time, (b) the subtraction of the fastest decay from the 
other curves, AA/ = M{tw = 7s) — A/(tu,),for ~ 
17.5(f), 25(T), 27(a), 35(»), 44(B). It can be seen that the 
M{tw=17) curve begins to overlap M{tw — 7) at approxi- 
mately 5000s. 



not a continuation of the initial stationary (power law) 
term. 

All measurements in this work were performed on an 
alloy CM.94Mn.06 that has a transition temperature Tg = 
31. 5i^. The sample has been the focus of man y pr evious 
studies and as such has been well-characterized l(f\ . Mea- 
surements were performed on a home built DC SQUID 
magnetometer that allowed for good control of the rapid 
temperature quench and subsequent waiting time. All 
measurements in this study were performed at T = 26K 
which corresponds to reduced temperature of 0.83Tg . 

We report very long time TRM decay measurements 
(10^ seconds) for very short effective waiting times ( 
7s < < 110s ). The data is plotted in Figure [TJi. The 
decays were obtained using identical fast cooling proto- 
cols with a small additional waiting time added before the 
magnetic field is shut off. It is this small additional wait- 
ing time that changes the age of the curves in Figure la. 
The ages of the individual curves were determined from 
the peak in the relaxation curve {S{t) — -~jjdM/dln{t)) 
plot shown in the inset of Figure It can be observed 
that as the effective waiting time increases the aging por- 



tion of the curves separate. This is a well-known effect of 
aging^3- We observe for the short waiting times curves 
that the decays overlap each other at very long times. 
Since all of the decay curves either follow or approach 
this long time decay, we conclude that it is independent 
of the waiting time and hence distinct from the aging 
regime. We plot the difference between the 7s waiting 
time curve and the longer waiting time curves in Figure 
^p. The range over which this difference goes to zero 
determines the onset of the third component of the de- 
cay for each of the waiting times. It can be observed 
that, within error, the decay of the 17s waiting time de- 
cay begins overlapping the 7s waiting time decay at ap- 
proximately 5,000 seconds. It can also be observed that 
the 25s waiting time decay overlaps, within error, at ap- 
proximately 35,000 seconds while the 27s waiting time 
decay overlaps at approximately 75,000 seconds. The 
35s and 45s waiting time decays systematically approach 
the 7s waiting time decay, but the magnetization dif- 
ference does not go to zero within a measuring time of 
100,000 seconds. This analysis is subject to the errors 
inherent within the measuring system. Sources of error 
include a high frequency point-to-point fluctuation and 
low- frequency "wobbles" (possibly due to small temper- 
ature fluctuations < IQmK). The signal to noise ratio, 
at 100,000 seconds, is approximately 100 to 1. 

It therefore appears that the aging effect is finite in 
time extent and at long enough time scale gives way 
to a long time component of the magnetization decay 
that is independent of the waiting time. Unfortunately 
the region of the overlap (5, 000-100, 000s) between the 7s 
and 17s waiting time decays gives only one and a half 
decades of measurement time with which to fit a func- 
tion, precluding an absolute determination of the form 
of the fitting function. We find that we can fit this range 
to a power law, a stretched exponential and a logarithmic 
function(Figure[21). Values of the fitting coefficients and 
exponents are given in the figure caption. We take a data 
point every 3 seconds for a total of approximately 20,000 
measurements in the overlap range. The large number 
of measured points and small signal to noise ratio give a 
very small value for the of the fitting functions. The 
stretched exponential is the worst fit with =5.8xl0~^. 
The power law fit gives =7.0x10"^, while the loga- 
rithm gives the best fit with =1.2x10^^, 

It is possible that Eq^Jis the correct description of the 
decay and that we are observing at long times, the short 
time stationary term after the aging decay has ended. To 
test this we have determined the upper limit of the power 
law exponent for the sample, measured at. 83 Tg, using a 
TRM magnetic field of 20G and a waiting time of 10,000 
seconds (data is displayed in[^). The short time decay in 
the TRM should be dominated by the power law decay 
of the short time stationary term, Eq^] Fitting over the 
first 3 seconds of this TRM decay we have determined 
that the upper limit for the exponent a is 0.01. This 
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compares with an exponent of the power law a = 0.072 
found for the very fong time decay. This value is also 
much larger than the value reported from ac susceptibil- 
ity experiments'! 7*1 for the short times stationary term. 
This significant difference suggests that the short time 
and long time decays are of different origin. 
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FIG. 2: Three fits to TRM(t„ = 7s), fitted to t > 5000s. 
The three fitting functions are: power law, stretched ex- 
ponential and logarithmic. For power law (M = Mot'"), 
Mo = 0.1881 and a = -0.072. For stretched exponen- 
tial (M = Moexp[f-°/ro], Afo = 1.48, tq = 10"^^ and 
a = -0.02559. For logarithmic fit( M = Mo + Log{t°')), 
Mo = 0.156 and a = -0.0148. 

We now address the question as to whether this long 
time decay is an additive term. Subtraction of any of the 
fitting functions from the general decay curves produce 
sets of curves that no longer scale as t/t^. In fact, we 
are not able to scale the subtracted data using any of the 
standard scaling techniques 0|. Unless we are to aban- 
don the concept that the aging decays as t/fm, or that the 
decays scale at all, we are led to the conclusion that this 
long time decay is not simply an additive term like the 
short time stationary decay. This analysis implies that 
this long time stationary decay is intrinsically related to 
the aging process and does not occur until the effects of 
aging are effectively over. An understanding of both, the 
mechanisms responsible for, and form of, this long-term 
decay can be found through analysis of the barrier model. 

The Barrier Model of the spin glass phase space was 
developed to explain aging in spin glasses. Inspired by 



the Parisi Til solution of the first principles Sherrington 
Kirkpatrick J_9J model of spin glasses and grounded in ex- 
perimental analysis, the barrier model implies that spin 
glass phase space can be described by diffusion of the 
system within a phase space, defined by a set of energy 
barriers which grow linearly as a function of Hamming 
distance|2l|. Generally a minimum barrier Aq is defined 
and the height of barrier N is equal to the product of N 
and Ao. The model simplifies the hierarchically ordered 
tree structure of the Parisi solution by collapsing local 
branches of the tree(states) into a single state with the 
appropriate degeneracy multiplication factor (this corre- 
sponds to the product of the barrier number N and the 
branching ratio r). Diffusion within the Barrier Model 
is essentially a trade-off between the geometrically grow- 
ing number of states as a function of increasing barriers 
and an exponentially decreasing probability for Boltzman 
hopping over the increasing size of the barriers. This 
model has proven to be very successful in quantitatively 
modelling aging decays, including temperature cycling. 
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FIG. 3: TRM experiment interpreted using a schematic of 
the Barrier model. During the cooling process, there is a 
"seeding" of the phase space. An initial distribution is set 
up. During the waiting time, states begin to diffuse into the 
phase space. The barrier associated with waiting time is Atw ■ 
The highest barrier is A^ax 

The initial state occupation is formed upon cooling the 
spin glass through its transition temperature to some 
measuring temperature, in a small magnetic field (Fig- 
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ure|31). Magnetic field cycling during the waiting time 
has previously been employed ,23j| to probe the redistri- 
bution of states in going from the field cooled manifold 
to the zero field manifold and back. The initial state oc- 
cupation was determined to be a delta function located 
at zero energy barrier with a possible background contri- 
bution from states close to the delta function. In reality, 
this study probed the initial distribution within the zero 
field manifold (after the field was cut) and this distribu- 
tion was also taken to be the initial distribution in the 
field cooled manifold after initial cooling. In the present 
study we find that it is not the case that the initial dis- 
tribution in the field cooled manifold is a delta function 
at the origin. Instead, we are only able to fit the long- 
time decay data by assuming that the states formed dur- 
ing the cooling process are distributed uniformly over all 
possible states (i.e. over all barriers). This has several 
significant implications. First, it has a self similar back- 
ground which is necessary for understanding the tem- 
perature dependence of the TRM. Secondly, it suggests 
that correlations of all possible strengths are formed dur- 
ing the cooling process. Finally, the uniform background 
produces a logarithmic decay, implying a limiting barrier 
height and hence a limit on aging. 

In Figure |21 we plot a fit of the whole 7s waiting 
time TRM decay using the results of the barrier model 
(Aq = .12, r = 1.048 and N^ax = 800 barriers), evolved 
from a constant background initial state, coupled with 
a power law term to model the "short time" stationary 
part of the decay (Eq. 1). Values for the power law 
were those determined from the 10,000 second waiting 
time data, previously mentioned. Extrapolating (straight 
line logarithmic dependence of the barrier model) to zero 
magnetization, we find a limiting age (corresponding to 
a maximum barrier) of t^ax ~l.lxlO^^ seconds. This 
corresponds to approximately 3000 years. This remark- 
able result implies that the phase space, at least in this 
sample, is limited in the time domain. 

We can test the assumption of an initial uniform state 
distribution over all barriers by performing an Isothermal 
Remnant Magnetization (IRM) measurement. In this ex- 
periment the sample was cooled in zero magnetic field, 
producing an initial distribution of states within the zero 
field (zero magnetization) manifold. The sample is kept 
in zero field for a time t^i — 10,000 seconds. A small 
magnetic field (20G) is then applied for a time tw2- After 
tw2 the magnetic field is removed and the decay observed. 
The data is plotted for several different values of m 
Figure ^ 

We start our analysis of the data with the assump- 
tion that the field cool (FC) and zero field (ZF) mani- 
folds are equivalent. This assumption has strong experi- 
mental backing. [23| The cooling procedure in the IRM 
should therefore produce an equivalent distribution of 
states in ZF manifold as the TRM produces in the FC 
manifold. Waiting for a long time t^i allows enough time 



for states, which contribute to aging, to diffuse over sig- 
nificantly large barriers in the ZF manifold. When the 
magnetic field is applied a small number of these aging 
states diffuse back to the zero barrier and jump to the 
FC manifold. These states then begin to age, starting 
from the zero energy barrier, in the FC manifold for a 
time tu,2- Removal of the magnetic field then produces a 
measurable decay of the states back from the FC man- 
ifold to the ZF manifold. If our assumption that the 
background is due to the cooling effect, then the uni- 
form background should be located deep within the ZF 
manifold and should not be observed in the IRM exper- 
iment. This is exactly what is observed. We do find an 
aging decay followed by a very small power law decay 
with an exponent compatible with that observed in the 
short times stationary measurements. 
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FIG. 4: IRM decays. We use t^i = 10000s and vary t™2 : 
10(0), 25(0), 100(A), 150(V), 300(0), 1000(<l) and 3000([>) 
seconds. The fitted curve is a power law with the power, 
a = -0.01. 

In conclusion, we have measured very long time TRM 
measurements with rapid cooling protocols and small 
waiting times. We find that that the aging time is fi- 
nite in extent and that there is a longtime decay which 
is independent of the waiting time. This decay is differ- 
ent from the short times stationary term and related to 
the mechanisms responsible for producing aging in spin 
glasses. 
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